Two-particle renormalizations in many-fermion perturbation theory: Importance of 

the Ward identity 
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We analyze two-particle renormalizations within many-fermion perturbation expansion. We show 
that present diagrammatic theories suffer from lack of a direct diagrammatic control over the phys- 
ical two-particle functions. To rectify this we introduce and prove a Ward identity enabling an 
explicit construction of the self-energy from a given two-particle irreducible vertex. Approximations 
constructed in this way are causal, obey conservation laws and offer an explicit diagrammatic control 
of singularities in dynamical two-particle functions. 

PACS numbers: 71.10Fd,71.10Li 



•4-* 



O 

(N 
> 

ON 

o 

(N 

o 

-I— > 



c 

o 
o 



X 



Correlated electrons in metals represent a system with 
unparalleled features that remain far from being satisfac- 
torily and fully understood. They are unique in that re- 
spect that to control their behavior in intermediate and 
strong coupling regimes we must have at our disposal ef- 
fective techniques for a reliable description of one- and 
two-particle characteristics. One-electron functions, such 
as the spectral function, the self-energy or the density of 
states carry information about the fermionic character of 
elementary excitations in metals. They are decisive for 
the Fermi liquid behavior as well as for deviations from 
the Landau quasiparticle picture. The two-particle cor- 
relation and Green functions control the critical behavior 
and signal cooperative phenomena and phase transitions 
via divergences in generalized susceptibilities. In inter- 
mediate coupling, where the Coulomb repulsion becomes 
comparable with the kinetic energy, dynamical fluctua- 
tions in the system become strong and may lead to ei- 
ther a phase transition to a new (magnetic) phase or to 
a breakdown of the Fermi-liquid behavior. In this transi- 
tion regime we have to treat both single-electron and pair 
excitations on the same footing and keep the two-particle 
functions directly accessible. 

One of the most flexible and physical ways to under- 
stand various phenomena of correlated electrons is to use 
many-body perturbation theory and Feynman diagrams. 
Except for the single-impurity Anderson and Kondo 
models, where the Fexmi-liquid regime survives to infinite 
interaction strengthjllu we have to renormalize the per- 
turbation expansion. Baym and Kadanoff were the first 
who showed how to implement renormalizations into t 
perturbation theory in a systematic and consistent way. 
The basic idea of their approach is to express physical 
quantities, and in particular the self-energy, as a func- 
tional of the renormalized one-particle propagator G and 
the bare interaction U, i. e., we construct a self-energy 
functional £ CT (k, iw n ) = S CT [G, U](k, ioj n ). Once we find 
an approximate form of this functional from a diagram- 
matic expansion free of self-energy insertions, we add the 
Dyson equation G~ 1 (k,iuj n ) = ^~ 1 (k,iw n ) - £ CT (k, iuj n ) 
to complete our approximation for the self-energy. We 
used Q to denote the bare one-particle propagator. If 



we extend our equations to situations with external per- 
turbing potentials, we can derive all necessary ther- 
modynamic functionso Thermodynamic consistence and 
macroscopic conservation laws are thereby guaranteed. 

Although we construct a functional for the self-energy 
from perturbation theory, the fundamental quantity in 
the Baym-Kadanoff approach is the generating functional 
related to the self-energy via a functional differential 
equation 



H a [G, U](k,iuj n ) 



5<S>[G, U] 
5G a (k,iu) n ) 



(1) 



When we are able to find the Luttingcr-Ward func- 
tional $[G, U] explicitly we speak about ^-derivable 
approximations.oEl 

Higher-order Green functions are derived from the 
self-energy functional via functional derivatives that may 
be viewed upon as generalized Ward identities. They 
connect lower-order with higher-order irreducible (ver- 
tex) functions. The identity connecting the two-particle 
irreducible vertex with the one-particle one (self-energy) 
reads in the direct space 



A^,(13,24) = 



<?S g (l,2) 

6 a G a ,(4,3) 6 a G a {2,l)5 a G a ,(4,3) 



(5 2 $[G, U] 



(2) 

where we denoted the space-time coordinates 1 = 
(Ri,ti) etc. We also introduced an index a denoting 
the appropriate two-particle irreducibility channel. We 
have three topologically nonequivalent two-particle ir- 
reducibility channels: electron-hole (a = eh), electron- 
electron (a — ee), and interaction (a = U) channels 
according to whether we cannot disconnect the diagram 
by cutting a pair of an electron and a hole, two electron 
(hole) propagators, or by cutting a polarization bubble, 
respectively, cf. Ref. |5} Using functional derivatives of 
the generating functional we derive all two-particle and 
higher-order Green functions. 

The formally exact approach of Baym and Kadanoff 
uses explicitly only mass renormalization, i. e., the self- 
energy functional depends explicitly on the renormalized 
one-electron propagator and the bare interaction. There 
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is no explicit two-particle renormalization in this formu- 
lation. The two-particle functions are passive outputs 
from functional derivatives of the self-energy functional. 
We hence cannot asses or directly control the critical be- 
havior via, for example, low-energy scalings or summa- 
tions of most divergent diagrams. Each change in the 
two-particle function arises only via an adequate change 
in the self-energy, which is too cumbersome and some- 
times even not viable. This is a severe drawback, in par- 
ticular in critical regions, where two-particle functions 
become singular and very sensitive to any small change. 

One can improve upon this by an explicit charge 
renormalization, that means that one replaces the bare 
Coulomb interaction in the perturbation expansion by 
appropriate two-particle irreducible vertices A". We then 
obtain a new generating functional in a form $[G, A Q ]. 
The most direct way to do this is to use the so-called 
parquet approach introduced in the nonrelativistic many- 
body theory by De Dominicis and Martini] The parquet 
approach differs from the Baym-Kadanoff construction, 
as standardly used, in that the parquet scheme takes the 
two-particle irreducible vertices as primary objects for 
which one tries to find a functional representation from 
the diagrammatic theory. 

Once we have prescriptions for the two-particle irre- 
ducible vertices A" we use the Bethe-Salpeter equation 
from the respective two-particle irreducibility channel to 
find the full two-particle vertex 

-{l + 5 aa ,)[h. a GGQ)T] aa ,{k;q,q'). (3) 

We introduced four-momenta k = (k, iu> n ),q = (q, iu m ) 
for fermionic and bosonic variables, with Matsubara fre- 
quencies w n = (2n + 1)ttT and v m = 2mirT at tem- 
perature T. Each two-particle scattering channel a is 
characterized by the way the irreducible vertices are in- 
terconnected by pairs of one-particle ^propagators GG, 
denoted here by the generic symbol 0fl 

Up to this point the one-particle propagator G and 
the vertices A Q have been treated as independent. To 
find a functional for the self-energy and then also for 
the generating functional $ one uses an exact Schwinger- 
Dyson equation of motion. Its explicit form depends on 
the model we choose. Here we use the lattice Hubbard 
model with completely screened, local Coulomb interac- 
tion. The Schwinger-Dyson equation then reads 

E CT (fc) = U^G-atf) - uY,r*-*(k; q, k' - k) 

k' k' q 

xG a (k + q)G- a (k' + q)G- <T (k'). (4) 

We introduced a short-hand notation = 

N ^ 1 J2q^ 1 J2m for 1 = (q^m) and analogously 
for k' — (k',iuv). The Schwinger-Dyson equation 
determines the self-energy functional T, a [G,A a ] from 
which we can construct a generating functional $[G, A"], 



which actually has been done for .various two-particle 
approximations of the parquet typeJj'l 

In the parquet approach all physical quantities are 
represented as functionals of renormalized one-particle 
propagators G a , representing mass renormalization, and 
two-particle irreducible vertices A^ CT , , standing for charge 
renormalization. The present way how the parquet-type 
approximations, or more generally two-particle renormal- 
izations, are applied, however, leads to a scheme effec- 
tively equivalent to that of Baym and Kadanoff. An in- 
tended direct control over the two-particle functions is 
lost at the end. We use equation (Q) to obtain the self- 
energy functional in the parquet approach. In conserving 
theories the physical two-particle functions must be de- 
termined via functional derivatives as in Eq. (^). The 
two-particle vertex constructed via a functional deriva- 
tive from the self-energy then differs from the vertex 
we started with. Furthermore, Eq. (Q) does no longer 
play the role of the Schwinger-Dyson equation, since the 
actual physical two-particle vertex T p h ys differs from T 
used in Eq.(^). This discrepancy reflects the general fact 
that we cannot fulfill both Eq. (||) and Eq. (0) unless 
we have an exact solution for the two-particle vertex. 
Hence, Eq. (Q) in the parquet construction serves as a 
generator of the self-energy functional from which all 
physical quantities are derived via functional differential 
equations. The two-particle vertices A a ,T we use in the 
parquet approach are only auxiliary functions and not 
the physical two-particle functions we need. 

It is clear that in approximate thermodynamically 
consistent and conserving schemes in the Baym-Kadanoff 
form we have to give up the Schwinger-Dyson equation of 
motion. A question arises whether also with two-particle 
renormalizations we have to lose the initial meaning and 
interpretation of the vertex functions. Actually, the two- 
particle renormalizations in the parquet approach were 
introduced in the effort to gain a better control over the 
behavior of two-particle functions in critical regions with 
singularities in Bethe-Salpeter equations. 

To regain the control over the two-particle functions, 
i. e., to return the meaning of two-particle functions to 
the diagrammatic representations for the vertices A Q , we 
have to resolve the Ward identity, Eq. (0), for the self- 
energy. The self-energy then will no longer be related 
to the two-particle vertex via Eq. (|j) but rather through 
an integral form of the Ward identity. In this way the 
two-particle functions will retain their complete diagram- 
matic representation enabling a direct control of their 
critical behavior. Using a Ward identity for the determi- 
nation of the self-energy from a given two-particle vertex 
guarantees thermodynamic consistence of the theory. 

The aim of this paper is to prove the following par- 
tially integrated Ward identity between the self-energy 
and the electron-hole irreducible triplet vertex 

z a (k)-z„(k') = J2K h Ak;q, k'-k) 

9 

x [G (J {k + q)-G u {k' + q)] . (5) 
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FIG. 1: Ward identity connecting the self-energy and the 
triplet irreducible vertex. The interaction acts vertically be- 
tween the upper and the lower line in two-particle functions. 
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FIG. 2: A typical self-energy diagram with internal loops (X) 
connected with the propagating electron only via the interac- 
tion (wavy) lines. 



Its diagrammatic representation is given in Fig. [j]. Equa- 
tion (|5|) holds for arbitrary four- momenta k and k' . It is 
a generalization of a Ward identity from noninteracting 
electrons in a random potential proved by Vollhardt and 
WolfleJa 

We stress that the Ward identity (||) can be proved 
only for the triplet irreducible vertex A CTCr . This is impor- 
tant in particular for the Hubbard model where the local 
Coulomb interaction acts only between particles with op- 
posite spins and singlet vertex functions A"_ CT are prefer- 
ably used. 

We use similar diagrammatic arguments as applied in 
noninteracting disordered systems to prove identity (||). 
Due to spin and charge conservation in the vertices of 
the perturbation theory, each diagram for the self-energy 
contains just a single electron trajectory connecting the 
incoming with the outgoing external electron line. We 



set aside this fundamental fermion trajectory propagat- 
ing the incoming charge and spin from the rest of the 
diagram that we denote X. Function X contains inter- 
action lines and only closed loops of fermion propagators. 
It is connected with the fundamental fermion trajectory 
via interactions as shown in Fig. |[ Each self-energy dia- 
gram can be classified according to the length (number of 
scattering events) of the fundamental fermion trajectory. 

We do not need to consider the Hartree term, since it 
does not contribute to the r.h.s. of Eq. (||). Hence the 
dynamical self-energy can be represented by an expansion 



£ CT (fc) = ]>>"+ 1 

n=l 

x ^ X^[G,U}(q 1 ,q 2 -qi,...,q„-q n -i,-q„) 



x G a {k + qi )G a {k + q 2 )... G a (k + q n ) (6) 

where n denotes the length of the fundamental fermion 
trajectory (n = 3 in Fig. |^). Note that the expansion on 
the r.h.s. of Eq. (^) is not an expansion in the interaction 
strength, since the loop function X contains interaction. 
It is clear that the loop function depends only on the 
transfer (bosonic) momenta q% , . . . , q n and not on the in- 
coming fermionic momentum k. 

Only the one-electron propagators contributing to the 
fundamental fermion trajectory are relevant for the proof 
of Eq. (|5|), since only there we can distinguish different 
external momenta on the l.h.s. of Eq. (||). We utilize the 
following identity 

G a {k + qi )...G a {k + q n ) - G a (k' + gi) . . . G a (k' + q n ) 

n 

= J2 G*(k + qi) ■ ■ ■ G a (k + g»_i) [G a {k + q t ) 
- G(j (k 1 + qi )} G a {k' + g i+1 ) . . . G (j {k! + q n ) (7) 
and rewrite the expansion for the self-energy to 



oo n 

£ CT (fc)-£ CT (fc') = ]Tt/" +1 ]r Yl X( n \q 1 ,q 2 -q 1 ,...,q n -q n _ 1 ,-q n ) 

ri— 1 i—1 qi , . ...q n 

x G(j {k + gi) . . . G(j {k + g i _i)G CT (fc' + q i+1 ) . . . G a (k' + g„) [G a (k + q { ) - G a (k' + qi)] 

OO n — 1 

= $>" +1 E A - - i - 1) (k;q,k'-k)[G <7 (k + q)-G (T (k' + q)] . (8) 

n=l i=0 

I 



We denoted contributions to an electron-hole vertex with 
i electron and n — i hole propagators from the fundamen- 
tal fermion lines of the two-particle function as A^™ ^ ■ 



The vertex A has the electron-hole structure, since the 
propagation in momentum k' is in the opposite direc- 
tion to the propagation of momentum k. The vertex A 
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must be irreducible in the electron-hole channel, since all 
reducible terms are already incorporated in the renor- 
malized one-electron propagators. On the other hand, 
each e/i-irreducible two-particle diagram can be uniquely 
closed to a self-energy via the fundamental fermion tra- 
jectory. The expansion in n and i hence covers all con- 
tributions_to the irreducible electron-hole vertex and we 
can writeE3 

oo n 

A£(k;q,q') = U £ U" +1 £ Afe"" (*; 1, </) • (9) 

n=0 i=0 

Inserting Eq. (||) into Eq. (g) we reveal identity (||). 

Relation (||) has two important consequences. First, 
continuity equation for two-particle correlation functions 
can be proved. Second, Eq. (|j) can be used to replace 
Eq. (Q) in the determination of the self-energy from a 
given two-particle vertex A^. 

To prove a two-particle continuity equation we use the 
following relation 

G a (k)G a (k + q) 

= G a (k)-G a (k + q) 

iv m - e(k + q) + e(k)-E a (fc + <?) + £, (fc) 1 ' 

in the Bethe-Salpeter equation for the two-particle 
function L aa (k;q,q') = G a (k)G a {k + q')[5(q) + 
rerer(fc; q, q')G a (k + q)G C r(k + q + q l ). The four-momentum 
delta function for q = (q, iv m ) reads 5(q) = NSq.oPSm^. 
Function L aa is the two-particle Green function from 
which the exchange term was removed so that it fulfills 
a Bethe-Salpeter equation. If we further utilize the sym- 
metry of the two-particle vertex we find another form of 
the Ward identity £ fe A q G a A e J l a {k; k'-k, q) = A q Z a (k'). 
When we multiply the Bethe-Salpeter equation for L aa in 
the electron-hole channel by the denominator of Eq. (|icj) 
and integrate over the incoming and outgoing fermionic 
momenta we obtain a continuity equation having in the 
limit of small transfer three-momenta q the following 
form 

iu m E aa (q, iv m ) - q • E£ CT (q, iv m ) = 0. (11) 



We denoted E rTrT (q,iv m ) = Y,k. g > L vcr{k;q' \q) and 
H^(q, iv m ) = J2 k ,q' Ve(k)L CTCT (fc; q', q). 

To use Eq. (|^) in the determination of the self-energy, 
we put k' = k but let the Matsubara frequencies iui n 
and iu) n i independent so that we can analytically con- 
tinue the self-energy difference to the case with u> + it] 
and ui — irj. We then use an analytically continued vertex 
function A^(k, oj + irj; q, C, 0, — 2ir/) in the Ward iden- 
tity ([s]) to construct the imaginary part of the self-energy 
S ff (k, to + irj) along the real axis. The intermediate com- 
plex frequency C takes values from an integration contour 
used to replace the sum over Matsubara frequencies. The 
shape of the integration contour depends on the analytic 
structure of the approximate vertex. The real part of the 
self-energy is calculated from the Kramers-Kronig rela- 
tion acomplishing thus a causal theory. This construction 
of the self-energy from the vertex function was already 
successfully applied in the parquet approach to noninter- 
acting disordered electron systemsJiil 

To conclude, the principal result of the paper is a 
partially integrated Ward identity, Eq. (||). It is a conse- 
quence of charge and spin conserving particle interaction 
and was proved by means of a diagrammatic expansion. 
The primary importance of identity (^|) lies in the pos- 
sibility to regain a direct diagrammatic control over the 
(critical) behavior of dynamical two-particle functions by 
defining the self-energy from the triplet vertex function 
via Eq. (||) . The present schemes with two-particle (ver- 
tex) renormalizations do not obey the Ward identity (|5|), 
except for its infinitesimal limit, and hence are unable to 
determine the self-energy from the physical two-particle 
vertices. Utilization of the Ward identity @ opens new 
possibilities to study critical behavior of correlated elec- 
trons, in particular, when we subject the electrons to a 
random potential and investigate the effects of electron 
interactions on the metal-insulator transition. 
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